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ABSTRACT 

Using a large sample of spiral galaxies for which 21 cm single-dish and/or long-slit optical spectra 
are available, we make a detailed comparison between various estimates of rotational widths. Differ- 
ent optical width estimators are considered and their limitations discussed, with emphasis on biases 
associated with rotation curve properties (shape and extent) and disk central surface brightness. The 
best match with Hi rotational velocities is obtained with Polyex widths, which are measured at the 
optical radius (encompassing a fixed fraction of the total light of the galaxy) from a model fit to the 
rotation curve. In contrast with Polyex widths, optical rotational velocities measured at 2.15 disk scale 
lengths deviate from Hi widths by an amount that correlates with the central surface brightness 
of the disk. This bias occurs because the rotation curves of galaxies are in general still rising at 2.15 
r<i, and the fraction of total mass contained within this radius decreases with increasing disk surface 
brightness. Statistical corrections, parameterized by the radial extent of the observed rotation curve, 
are provided to reduce Polyex and Hi width measurements into a homogeneous system. This yields a 
single robust estimate of rotational velocity to be used for applications of disk scaling relations. 
Subject headings: galaxies: kinematics and dynamics — galaxies: distances and redshifts — galaxies: 
spiral — galaxies: structure — cosmology: observations 



1. INTRODUCTION 

The importance of rotational velocities of spiral galax- 
ies for extragalactic astronomy and observational cosmol- 
ogy s tems in large part fro m their use in the Tully-Fishcr 
(TF; iTullv fe Fisherlll977f ) relation. TF surveys initially 
relied on rotational velocities measured from global (i.e., 
spatially integrated) Hi line profiles, typically obtained 
using the Arecibo, Nancay, Green Bank, Parkes, and Ef- 
felsberg radiotclcscopes. Owing to sky coverage or sensi- 
tivity limitations of single-dish antennas, optical long-slit 
spectroscopy has been widely used to complement 21 cm 
Hi observations for TF applications. The accuracy of 
the velocity width measurement is particularly impor- 
tant for a TF survey, since the associated uncertainties 
play a do minant role in the overa ll error budget of the 
technique (jGiovanelli et al.lll997bl hereafter G97b). 

The observed velocity width of a spiral galaxy spec- 
trum is primarily determined by the shape of its rota- 
tion curve (RC), the distribution of the adopted kine- 
matic tracer, the inclination of the plane of the disk, 
and the measurement method. Integrated, single-dish 
Hi profiles provide very limited information regarding 
the distribution of the Hi gas, the shape of the RC, the 
presence of disk asymmetries, departures from planarity 
or non-circular motions. However, the Hi gas is largely 
distributed in the outer parts of disks, sampling pref- 
erentially regions where the RC is flat or slowly rising; 
in addition, these measurements are completely indepen- 
dent of uncertainties in the determination of the disks' 
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projected major axis. Long-slit, major axis Ha spectra 
yield an estimate of the shape of the galaxy's RC, and 
they do so with substantial angular resolution, allowing 
one to gauge departures from circular motion and disk 
asymmetries. Furthermore, they can probe the inter- 
nal kinematics of galaxies in high density environments, 
where the Hi signal could be confused by the presence 
of more than one target within the telescope beam, or 
where the Hi emission could be truncated because of in- 
teractio n with the intracluster medium (Hi deficiency ; 
see, e.g.. lHavnes fc Giovanellil[i984t ISolanes et alJ .2001). 
Long-slit optical RCs are also easier and faster to ob- 
tain at higher redshifts (i.e., cz larger than a few thou- 
sand km s -1 ) compared to Hi spectra. However, they 
only sample the actively star forming regions of the disk, 
which typically extend only about half as far out as the 
detectable Hi gas, and are sensitive to errors in the es- 
timate of the major axis position angle, high internal 
extinction, and other resolution-related effects (these is- 
sues are particula rly important for int ermediate redshift 
galaxies; see, e.g.. lKapferer et aT1l2006| ). Thus widths ob- 
tained with the two methods are differently affected by 
systematic effects. 

The cross-calibration of Ha and Hi velocity width 
scales and the understanding of the systematics involved 
in each technique are essential for the combined use of 
heterogeneous data sets in TF work. In fact, system- 
atic differences between measures of rotational velocities 
derived separately from optical and radio spectra could 
mimic spurious bulk motions in peculiar velocity surveys, 
if the samples are segregated in the sky or in redshift 
space. Such width conversion is also important for kine- 
matic studies of intermediate redshift (z > 0.4) galax- 
ies that must rely on other tracers, such as [O II]A3727. 
Widths measured with these diagnostics must be cali- 
brated with tho se obtained with Hi and Ha spe ctroscopy 
of local spirals (jKobulnickv fc Gebhardll2000l ). 
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As discussed in the first paper of this series 
(jCatinella. Haynes. fc G iovanelli 2005|, hereafter Paper 
I), the choice of the best algorithm to measure rota- 
tional velocities from long-slit spectra is not an obvious 
one. Therefore here we compare optical widths obtained 
from velocity histograms and from model fits to the RCs 
and, for the latter, address the issue of the choice of 
the reference radius for width measurement. In par- 
ticular, we discuss the vulnerability of optical widths 
to systematic effects associated with RC shape and ex- 
tent and disk central surface brightness (SB). For the 
RC width estimator that is least affected by these bi- 
ases, we provide the cross-calibration with single-dish 
Hi widths. For a discussion of how rotational veloc- 
ities derived from global Hi profiles may differ from 
thos e derived fr om H i RCs, the reader is referred to, 
e.g., IVerheiienl (|2001h . Comparisons between Hi and 
Ha velocity widths have been previously carried out by 
other a uthors, including Mathcwson . Ford, fc Buchhornl 
(|1992t l. [Ravchaudhurv et al.l (|1997f ). and lCourteaul (|1997l . 
hereafter C97), based on overlap samples substantially 
smaller than the one used here. However, the conversion 
relations suggested by previous authors do not take into 
account the biases discussed in this work. In particular, 
we show that the cross-calibration of optical and radio 
widths depends on the shape and extent of the observed 
RCs and, for widths measured at 2.15 disk scale lengths 
(e.g., C97), on the SB of the galaxies. 

In Paper I we presented the latest installment of 
Ha s pectroscopy for the SFI+- (- sample (jSpringob et al.1 
2007) and described our technique to measure rotational 
velocities from long-slit spectra. We also showed that 
RC extent and Hi content are correlated, in the sense 
that for a given optical size, galaxies with larger Hi con- 
tent are characterized by more extended RCs. This re- 
sult, which holds for objects with normal Hi content (i.e., 
not only for Hi-deficient, cluster members), will be used 
here to explain observed trends in the comparison be- 
tween optical and radio widths. The average, or tem- 
plate, RCs of disk galaxies binned in luminosity int ervals 
obtained bv lCatinella. Giovanelli. fc Haynesl (|2006l here- 
after CGH06), together with the SB effects on observed 
disk scale lengths and optical sizes discussed therein, are 
also relevant for the analysis presented here. 

The plan of this paper is as follows. The sample and 
velocity width measurements from optical RCs and Hi- 
line profiles used in this work are described in S}2j The 
comparison between optical and radio rotational veloc- 
ities is presented in Ej3j Systematic and SB effects are 
discussed in $4] and our conclusions summarized in £0 

2. DATA SAMPLE AND VELOCITY WIDTH 
MEASUREMENT 

Our data set, referred to as SFI++, is a homoge- 
neous compilation of /-band photometry and long-slit 
Ha and/or Hi rotational parameters for ^5000 galaxies, 
mostly late spirals. The SFI++ catalog was designed for 
peculiar velocity studies in the l ocal universe us i ng th e 
TF method and is presented in ISpringob et al.1 ([2007) . 
Most of the SFI++ galaxies (85%) have velocity in the 
CMB frame cz<10,000 km s -1 (97% have cz<15,000 
km s -1 ). Corrected velocity widths range between 53 
and 834 km s -1 (and between 200 and 500 km s -1 for 
80% of the objects), with a mean of 331 km s — 1 . 



The sample used in this work consists of 2753 SFI++ 
galaxies with high quality long-slit Ha/[N II] RCs and I- 
band photometry, for which good quality Polyex fits (see 
£|2.ip are also available. Disk central SB values, (1q (mea- 
sured from the /-band profiles and corrected for Galactic 
and internal extinction, cosmological k-term, and con- 
verted to face-on perspective), for this sample vary be- 
tween 14.0 and 21.3 mag arcsec -2 (only two objects have 
/xq < 14.0 mag arcsec -2 ); the distribution of ^ peaks 
at 18.6 mag arcsec -2 . This data set was used for the 
comparisons involving optical widths only. For those in- 
volving radio widths, we used the subset of galaxies with 
high quality Hi spectra (uncertain or confused Hi profiles 
were discarded, as well as a few outliers whose optical 
and radio widths are clearly inconsistent, yielding a final 
sample of 873 objects). 

2.1. Optical Widths from Long-slit Rotation Curves 

The algorithms commonly used to measure rotational 
widths from long-slit spectra are based on either velocity 
histograms or functional fits to the RCs. For the for- 
mer, the width is obtained as the difference between two 
percentile points of the RC velocity histogram; for the 
latter, the width is measured at a fixed radial distance in 
the disk (obtained from photometry) from the RC model 
fit. The histogram technique is simple, resembles that 
used to measure rotational widths from Hi integrated 
line profiles (although not intensity weighted), and does 
not require accompanying photometry. However, all the 
spatial information contained in the RC is lost and sys- 
tematic biases arise as a consequence. In this work we 
characterize the biases affecting histogram width deter- 
minations and compare model widths measured at two 
different radial distances. We adopt the following RC 
width definitions: 

(a) Wso — Vgo — Vio, the difference between the 90th 
and 10th percen tile points of the velocity histogram 
(e.g.. IVogtl |1995[ ) . We determined Vio and Vgo for our 
RCs by linear interpolation between the nearest velocity 
data points. This definition is preferable over the dif- 
ference between th e maximum and the minimum veloc- 
ity of the RC (e.g.,lMathewson. Ford, fc Buchhornlll992t 
iMathewson fc F ord 1996|), because it is less sensitive to 
the noise level of the RC (the maximum and the min- 
imum velocities typically correspond to the outermost 
detectable Hn regions, where the Ha emission is usually 
fainter). 

(b) Wpe- This is our algorithm of choice for RC width 
measurement and is described in detail in Paper I. In 
summary, we fold each RC and fit the following function 
(Polyex model) to it: 

VpeM = Vb(l - e-^Xl + ar/r PE ) (1) 

where Vo, tpe, and a determine the amplitude, the ex- 
ponential scale of the inner region, and the slope of 
the outer part of the RC, respectively. The Polyex 
width is measured from the fit at the optical radius 3 

3 Ropt is defined as the radius encompassing 83% of the total 
integrated light. It is equivalent to 3.2 disk scale lengths for a Free- 
man disk, i.e. an exponential disk with central surface brightness 
HO = 21.65 mag arcsec -2 in B-band (corresponding to fia ~ 20 
mag arc sec -2 in /-band , adopting an average B — I color index of 
1.7 from lde JongllT996T) . 
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R op t' W / pe=2 VpE(-Ropt)- Other fitting functions have 
been proposed, either e mpirical (e.g., C97; lKratsov et all 
119981 IVogt et all I2004D or based on physical m odels of 
disk galaxies (e.g.. iPersic. Salucci fc Stel fl996h . The 
Polycx model has no physical motivation, but has enough 
flexibility to fit the vast majority (~94%) of the individ- 
ual RCs in our data set, as well as the average RCs in 
luminosity intervals presented in CGH06. For a < 0, it 
also allows one to fit declining RCs. 
(c) Wz.2- We also discuss model widths measured from 
the Polyex fit at 2.15 (hence W2.2) disk scale lengths r^, 
another common choice of reference radius for RC width 
measurement (see §2.1. 1|) . 

Observed velocity widths derived from RCs must be 
corrected for cosmological broadening (to obtain the rest- 
frame velocities) and deprojected to an edge-on view: 



where z and i are redshift and inclination to the line of 
sight of the galaxy, respectively. The uncertainties in- 
troduced by the inclina tion correction are d iscussed in 
IGiovanelli et al.l (|1997al hereafter G97a) and iDale et all 
(|1997f h Corrections for turbulent (i.e., non-circular) mo- 
tions of the gas, usually applied to widths derived from 
Hi-line spectra (see H2.2§ . are considered negligible for 
RC widths. The reason is that, at any spatial location 
along the slit, the RC is obtained from the velocity of 
the peak of the line emission, and the broadening of the 
line due to turbulent motions only affects its width (more 
details on RC extraction can be found in Paper I). Two 
other effects that could distort the observed RCs are slit 
smearing (because of the finite width of the spectrograph 
slit, regions that are slightly off the major axis could con- 
tribute to the Ha flux, smearing the RC) and slit mis- 
alignment. These effects preferentially act in the sense 
of underestimating the true velocity widths. Numeri- 
cal simulations (to be presented elsewhere) show that 
the combined effect of slit smearing and misalignment 
has a typical amplitude of the order of 1-2%, signifi- 
cantly smaller than the accuracy of width measurements. 
Therefore we do not correct the observed widths for these 
effects. 

Velocity widths used in this work are corrected for cos- 
mological broadening only; the sin i correction is not 
applied in order to minimize the uncertainties in our 
analysis, but we checked our results for possible depen- 
dences on inclination. The plots shown in this paper do 
not change significantly when the inclination-corrected 
widths are used instead. 

2.1.1. Which Reference Radius for RC Width Measurement? 

The choice of the radius at which the model velocity 
widths should be measured is subject to some debate. As 
discussed in Paper I, observed RCs usually do not show 
a feature that could be easily associated with the maxi- 
mum circular velocity. In particular, RCs are in general 
flat or even rising beyond 2.15 ra, i.e. the lo cation of 
the p eak velocity for a pure exponential disk (|Freemanl 
1970). Compared to rotational velocities measured at 
2.15 ra, widths evaluated at i? pt are less affected by ex- 
tinction and are more resilient, in the sense that a given 
error on the spatial scale translates into a smaller scatter 



in the measured widths, because velocity gradients are 
smaller in the outer disk. Moreover, i? pt is a cumula- 
tive parameter, and hence more straightforward to mea- 
sure than rd and less subject to extinction effects. The 
use of i? pt avoids an additional complication. When 
the reference radius is a photometric quantity, such as a 
multiple of or an isophotal radius, its location on the 
disk will vary with the inclination of the galaxy to the 
line of sight. As a result, the corresponding measured 
width will in general depend on inclination in a more 
complex form than purely sinusoidal. The process could 
be standardized by using the radius at a fixed inclination, 
e.g. face-on. However, the measured scale length of 
a given galaxy varies by about 50% purely as the disk 
incli nation varies from int ermediate to high values (see 
e.g. IGiovanelli et al.lll994 Fig. 12), and the conversion 
r c i — » r^j is quite uncertain. Thus, the reference radius 
should be chosen in the very outer disk, where one may 
safely assume nearly total transparency and the conver- 
sion to a face-on value is straightforward. 

One advantage of width measurements at 2.15 is 
that most observed RCs are mapped as far out as this ra- 
dial distance, whereas only ~60% reach i? op t- Thus, the 
measurement of Wpe at i? pt requires more frequently 
an extrapolation of the fit. However, such extrapolation 
is mild, because, even if the RCs still rise near i? pt, they 
do so at a small rate. Moreover, most of the RCs in our 
sample (84%) extend to at least 0.8 i? op t (only 2% have 
extent smaller than 0.5 i? pt- These are not used for TF 
applications). 

For all the above reasons and in our experience, the to- 
tal uncertainty associated with the width measurement 
at i? pt is smaller than that at 2.15 r^- Ultimately, opti- 
cal widths need to be compared to Hi ones, which sample 
the galaxy kinematics out to larger radial distances in the 
disks. These should therefore provide more reliable es- 
timates of the maximum rotational velocity, at least for 
objects in low density environments, where the Hi distri- 
bution is undisturbed. Such a comparison is presented 
in §0 

2.2. Radio Widths from Hl-line Profiles 

Observations with single dish radio telescopes, whose 
beam areas significantly exceed the angular size of the 
Hi distribution of a target galaxy, yield a spectrum of 
the galaxy's spatially integrated Hi profile. Such profiles 
reveal neither the shape of the RC nor the details of the 
Hi distribution, an argument that has sometimes been 
used to suggest limited usefulness for single dish 21cm 
widths (e.g. C97, §11). Nonetheless, such widths yield 
TF relations of low scatter. This occurs because of the 
small slopes of RCs in the outer parts of disks, where 
much of the gas tends to be located (typically half of 
the Hi lies outside of about three scale lengths from the 
center). This produces sharp outer edges to the spectral 
profiles, so that a width measurement can be obtained 
with high precision (generally <10 km s" 1 ). This cir- 
cumstance is more significant for luminous objects. For 
galaxies with rotational velocities of ~ 75 km s _1 or less, 
the outer slope of the RC is large and the turbulent mo- 
tions become relatively more important. As discussed 
in G97b, the TF scatter among less luminous objects is 
quite large, partly because the shape of the RC and the 
presence of turbulence make the derivation of rotational 



4 



Catinella, Haynes, & Giovanelli 



velocity more uncertain, as are estimates of disk inclina- 
tion, but principally bec ause the TF intrinsic scatter is 
larger among them (see iHoffman et al.l [1996 as well as 
G97b). Therefore, these objects are ill-suited as targets 
for TF studies for the determination of distances or pecu- 
liar velocities, and such studies may be better off without 
them. 

Rotational widths of disk galaxies are commonly ob- 
tained from double-horned Hi profiles by measuring 
the velocity difference between the outer N% levels of 
each peak. Our techn ique is described in detail in 
lHavnes et all (|1999f > and lSpringob et al.l (|2005l hereafter 
S05). In summary, a straight line is fitted to the receding 
side of the Hi profile, between the 15% and 85% levels of 
the peak flux; V r is the velocity for which the line fit has 
a flux equal to 50% of that of the peak. The analogous 
process on the approaching side of the profile yields V a ; 
thus the measured width is W / 2i= V T — V a . 

Standard corrections applied to Hi-line widths include 
those for instrumental and cosmological broadening, tur- 
bulent motions, and deprojection to edge-on view: 

/ ^Obs _ A \ J 

W co " = — - A t — (3) 

\ 1 + z J sin i 

The cosmological and inclination corrections are the 
same as those applied to the optical widths (see eq. [2]). 
The term A s accounts for the effects of the spectrome- 
ter resolution, the amount of smoothing applied, and the 
signal-to-noise quality of the spectrum; the exact func- 
tional dependence of A s on the latter two quantities must 
be modeled with numerical simulations (S05). This cor- 
rection is of the order of the width of the velocity chan- 
nels, typically between 5 and 11 km s _1 for our data. 
The correction for turbulent motions, A t , is somewhat 
uncertain, and different values have been adopted in the 
literature. S05 derived it for our data by introducing 
a random, isotropic o =10 km s^ 1 turbulent motion to 
their simulations, and obtained A t =6.5 km s _1 , inde- 
pendent of profile width and signal-to-noise ratio. 

As for the optical data, no inclination correction is 
applied to the Hi widths. 

For convenience, the definitions of all the radial 
distances and rotational velocities adopted in this work 
are summarized in Table [TJ 

3. VELOCITY WIDTH COMPARISONS 

3.1. Optical Widths 

As mentioned in £12.11 width determinations based on 
velocity histograms do not make use of the spatial infor- 
mation contained in the RCs and this causes systematic 
deviations with respect to measurements performed at a 
fixed radial scale. On average, since the RCs are typi- 
cally rising within the optical disks (CGH06), histogram 
widths are expected to be systematically smaller (larger) 
than Wpe for RCs less (more) extended than i? ptj where 
the Polyex width is evaluated. This is shown in Figure 
[TJ where we plot the difference between Polyex and his- 
togram widths as a function of RC extent (expressed in 
units of i? pt ) ■ The effect is smaller for RCs with larger 
extent, because galaxies with intrinsically large Ha disks 
are also more luminous, and therefore their RCs tend to 
be flatter (CGH06). 



In order to characterize the dependence of histogram 
width measurements on the velocity gradient of the RCs, 
we divided our sample into four bins with increasing RC 
outer slope, as obtained from the Polyex fit between 0.5 
and 1.0 -R pt- Figure^ shows the difference between Wpe 
and Wsq as a function of RC extent for the four bins, 
roughly corresponding to falling (a), flat (b), mildly ris- 
ing (c), and rising (d) RCs. Several features can be no- 
ticed: (1) the difference Wpe — Wgo depends on both 
extent and slope of the RC; (2) Wso and Wpe yield the 
same results only for flat RCs, traced out to a radial dis- 
tance of the order of i? pt; (3) over the -R ma x/-Ropt range 
considered here, W$q is smaller than Wpe for rising RCs, 
the difference increasing for RCs with steeper outer slope 
and smaller extent. In the case of falling RCs, Wgo is 
larger than Wpe; (4) the intercept of the relation with 
the axis -R m ax/-Ropt = l increases from approximately —20 
to 20 km s _1 from falling to rising RCs. 

We performed numerical simulations to model the 
trends observed in Figured) We simulated mock catalogs 
of galaxies with ideal, perfectly symmetric and noiseless 
RCs, assigning RC extent (R m ax/R op t) and Polyex fit- 
ting parameters with the same frequency distributions 
observed in the SFI++ sample (paying attention to the 
fact that the distributions of the RC parameters are not 
independent from each other). We calculated RC slopes 
and widths, and binned the simulated data as done for 
our observations. The results are shown as solid lines 
in Figure [2] Even with the simplistic RC description 
adopted, our simulations are able to reproduce the sys- 
tematic trends observed in the data. A more realistic 
modeling, taking into account statistical distributions of 
Ha emission extent, RC asymmetry, and measurement 
errors, could easily explain the remaining (small) differ- 
ences between the data and the results of the simulations 
presented here. These results demonstrate that the sys- 
tematic biases affecting histogram width measurements 
are completely explained, once the extents and shapes 
of the RCs are taken into account. RC widths based on 
velocity histograms will not be further discussed. 

The next issue we want to address is the choice of 
the reference radius at which model widths should be 
evaluated, i? pt or 2.15 r<j. The uncertainties affecting 
the measurement of i? pt and r<j have been discussed in 
H2.1.11 in this section and in H3.2I we will compare Wpe 
and W2.2 with each other and with velocity widths de- 
rived from Hi spectra. 

The mean width difference Wpe — W2.2 for our sample 
is 17 km s _1 (with a dispersion er=19 km s^ 1 ). This 
offset is expected from the fact that, on average, the 
RCs are still rising at 2.15 r d (CGH06, Fig. 1) and W PE 
is measured farther out in the disk. However, values of 
Wpe — W2.2 for individual galaxies show a large variation, 
from —32 to over 100 km s _1 (see Fig. Obelow). 

As opposed to histogram widths, model widths do not 
depend on RC extent, as long as the RC is sampled to 
a radius large enough to warrant a reliable fit. In fact, 
Wpe — W2.2 does not show systematic variations when 
plotted as a function of Rmax/Ropt (not shown). On the 
other hand, the difference between widths measured at 
-R op t and 2.15 rd should increase for larger i? op tAd ratios, 
and thus should depend on SB, because galaxies with 
larger values of Ropt/fd are typically of higher optical 
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SB (CGH06). 

In order to quantify these biases and gain more insights 
on the role of the SB of a galaxy, we studied the varia- 
tion of the width difference as a function of the i? pt / fd 
and -R23.5/fd ratios 4 and of the /-band disk central SB, 
/XQ. Figure (3ja, shows the dependence of Wpe — W2.2 on 
the i?opt/^d ratio for the whole sample. This correla- 
tion is entirely expected from the fact that the RCs are 
typically still rising at 2.15 r d , where W2.2 is measured. 
Therefore, for a given R opt , Wpe — W2.2 will increase 
for smaller values of r<j, i.e. larger -RoptAd ratios, as 
observed. The bins corresponding to larger width differ- 
ences have larger error bars, due to both smaller number 
statistics and larger scatter among the individual data 
that were combined. 

As already mentioned, galaxies of higher /i have typ- 
ically larger values of i? op t/ r d and thus, based on Figure 
[3k., they should also have larger Wpe — W2.2 width differ- 
ences. This is confirmed in panel (b), where the sample 
has been divided into three intervals of central SB (cor- 
responding to those adopted in Fig. 9 of CGH06, except 
that the two central bins in that figure have now been 
combined). 

It is perhaps more insightful to study the dependence 
of the width difference on the R23.5/r d ratio, since an 
isophotal radius has a more straightforward dependence 
on the disk central SB than i? pt- In fact, the observed 
brightness profile of an exponential disk, expressed in 
magnitude terms, can be written as: 

(i(r) = no + 1.086 r/r d 

which implies a linear relationship between no and the 
#23.5/^ ratio: 

^0 = 23.5 - 1.086 R 23 . 5 /r d 

Since i? pt and -R23.5 are clearly correlated, this also im- 
plies a correlation between no and Rapt /fa, but the scat- 
ter is of course larger. To this extent, it is worth remem- 
bering that -R pt is defined as the radius encompassing 
a fixed fraction of the total light of the galaxy, i.e. it is 
not a pure disk quantity. Therefore, its value depends on 
the degree of concentration of the light as well as on the 
possible presence of a spheroidal component. For a given 
total luminosity, smaller exponential disks clearly have 
larger i? ptAd ratios. For a Freeman disk, R op t/r d = 3.2 
(see footnote 1); galaxies of higher SB and similar val- 
ues of r d have larger i? op t/^d ratios. However, notice 
that a galaxy harboring a bulge has a smaller R op t/r d 
than a galaxy with the same exponential disk but no 
bulge. Figures [3b and [3fcl illustrate the behavior of the 
width difference as a function of i?23.sAd for the whole 
sample and separately for the three SB intervals, respec- 
tively. Panel (d) shows that galaxies of different SB are 
well segregated in the i?23.sAd parameter space and that 
larger width differences are clearly associated with higher 
disk central SB. This will be further discussed in £j4j 

In order to decide which model width, Wpe or W2.2, 
is best suited for use in the TF relation, we will com- 
pare these with width measurements obtained from Hi 

4 ^?23.5 is the isophotal radius measured at an /-band surface 
brightness of 23.5 mag arcsec -2 . This approximately corresponds 
to R2 5 in B-band, a ssuming an average B — I color index of 1.7 
ffrom ldTjoiTari99^ . 



spectra in ^3.21 The optical width of choice will be that 
providing the best match to W21 based on the amount 
of offset and scatter of the correlation and, most impor- 
tantly, vulnerability to systematic effects. 

3.2. Optical vs. Radio Widths 

The results presented in the previous section showed 
that the difference between optical model widths mea- 
sured at i? op t and 2.15 r d depends on the SB of the 
galaxy. Thus we first evaluate the impact of SB-related 
biases on comparisons between radio and optical rota- 
tional velocities by plotting the corresponding width dif- 
ferences as functions of the i?23.s/^d ratio. As Figured] 
clearly demonstrates, Wpe provides a more reliable esti- 
mate of the full width of disk galaxies, as measured from 
spatially integrated Hi-line profiles, than does W2.2- In 
fact, not only is the difference between W21 and Wpe 
always smaller than that between W21 and W2.2, but 
it is also not affected by systematic trends, whereas 
W21 — W2.2 increases for larger i?23.5Ad ratios, i.e. galax- 
ies of higher SB. Plotting the width differences as func- 
tions of the central SB yields analogous results. There- 
fore we will only cross-calibrate Hi and Polyex widths 
and discuss the implications of the SB effects associated 
with RC widths measured at 2.15 r d in gj 

As seen in Figure[4] Hi widths are systematically larger 
than Ha ones; the average value of W21 — Wpe is 14 ± 1 
km s _1 (with a dispersion a=29 km s _1 ). Interestingly, 
the width difference tends to increase for galaxies with 
extended RCs, as shown in Figure [5k- The velocity in- 
crease is small for the whole sample and the scatter is 
large (excluding the two bins with -R ma x/-Ropt>2, which 
include less than 20 objects each, the velocity dispersion 
increases monotonically from 20 to 53 km s _1 with de- 
creasing RC extent), nonetheless the effect is systematic. 
The same bias can be seen in panel (b), where the ratio, 
rather than the difference, of the two widths is plotted. 
The solid lines in panels (a) and (b) are quadratic fits to 
the data; their analytical expressions can be found in the 
first two rows of Table [U respectively. 

The fact that the RCs might still be rising beyond i? op t 
could explain the offset between W21 and Wpe (as long 
as Hi emission is detected beyond i? op t), but not the 
dependence of their ratio on RC extent. The latter re- 
quires that RC and Hi extents are also correlated. In 
Paper I we showed that, for a given optical size, star for- 
mation can be traced farther out in disks of galaxies with 
larger Hi content. We also pointed out that this result 
is not limited to the well known case of truncated star 
formation in Hi-poor cluster objects, but also holds for 
Hi-normal, field spirals. From the analysis of a sample 
of 108 Hi-rich galaxies with available Hi synthesis maps, 
iBroeils fc Rheel (|1997f) showed that Hi mass and Hi di- 
ameter are tightly correlated; moreover the Hi diameter 
correlates well with the optical size of the galaxy. There- 
fore, a correlation between RC extent and Hi size seems 
to be a reasonable assumption. Thus, biases might arise 
because for galaxies with larger Ha extent the Hi width 
traces farther out into the halo, while we continue to mea- 
sure Wpe at i? pt ■ This also implies that the dependence 
of the width ratio (or difference) on RC extent should be 
larger for rising RCs and negligible for flat ones. This 
is clearly demonstrated in Figure O where the sample 
is divided into two classes with different RC slopes at 
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the optical radius. The dependence of W^i/Wpe on RC 
extent is present for galaxies with rising RCs (a), but 
effectively disappears for the rest of the sample (b). Lin- 
ear fits to the data points are shown as solid lines and 
listed in the last two rows of Table O Notice that, for 
rising RCs with extent smaller than i? op t, Wpe tends to 
be larger than W21, implying that the extrapolation of 
the RC to i? op t (where Wpe is measured) often overes- 
timates the width measured from Hi profiles. However, 
it should also be noticed that the first bin of panel (a) 
includes only 20 elements and that the scatter in these 
plots is large (the value of the la dispersion around the 
mean is 0.15 for the first bin and ~0.10 for the other bins 
of panel [a], and varies between 0.07 and 0.20 for panel 
W). 

The common assumption that optical widths underes- 
timate radio ones, because RCs are often rising and the 
Hi gas extends farther out in the disk compared to Ha 
emission, is confirmed by our data. As argued above, 
the observed dependence of W21/WPE on RC extent is 
also in agreement with this picture, if one assumes that 
galaxies with larger RC extent also have larger Hi size. 
However, these arguments do not explain why Hi widths 
are systematically larger than Ha ones for galaxies with 
flat RCs (Fig. [BJa). Even adopting a much stricter defi- 
nition of flat RC (i.e., —0.1 < < 0.1, a requirement 
met by 161 objects) yields an average velocity offset 
W 2 i - Wp E of 17 ± 2 km s^ 1 . At least part of this 
offset might indicate that the turbulent motion correc- 
tion applied to radio widths (6.5 km s _1 , see § 12. 2|) is 
underestimated. Not only is such a correction somewhat 
uncertain, with various values being advocated by differ- 
ent authors, but also it does not take into account tidal 
disturbances, bars, and non-planar motions which can 
occur frequently and with varying amplitude. Another 
important factor to keep in mind is that widths obtained 
from global Hi profiles depend on both RC shape and 
distribution of the Hi gas. Resolved Hi RCs of nearby 
galaxies show that, for most intermediate and bright ob- 
jects, a velocity peak is often reached before or near 
~3 t a , followed by a small decline and a flat outer part 
(e.g-ICasertano fc van GorkonJ[l99lt lBroeilslll992t C97; 
IVerheijenll2001l ). In these cases, the width of the global 
Hi profile measures the peak velocity, whereas the Polyex 
width obtained from the Ha RC would yield the ampli- 
tude of the flat part (as typically i? op t/rd>3.2 for high 
SB galaxies). 

The results presented in this section indicate that, in 
order to cross-calibrate optical and radio rotational ve- 
locities, both extent and shape of the RCs should be 
taken into account. The relations presented in Table [2] 
allow us to combine Polyex and Hi velocity widths into 
a homogeneous data set, which can be used for appli- 
cations of the TF distance method without introducing 
significant systematic biases. 

4. DISCUSSION 

4.1. Systematic Effects 

The comparison between Hi and Ha veloci ty widths 
has b ee n carried out by several autho r s. | Courteau 
(119921) , iMathewson. Ford, fc Btichhornl (|1992f L IVogi 
(1995), and iRavchaudhurv et al.l (|1997t ) found tight cor- 
relations between radio and optical histogram widths, 
with slopes approximately equal to unity, small zero 



points (< 5 km s _1 ), and dispersions ranging between 
10 and 30 km s . The most comprehensive calibration 
to date is that of C97, who compared optical widths mea- 
sured with several algorithms for a sample of 304 galaxies 
with Hi widths obtained by G97a (we note that the G97a 
data set is part of the SFI++ sample, but the Hi widths 
have recently been recomputed using different instrumen- 
tal and turbulent motion corrections, following S05). He 
concluded that the best match to 21 cm widths is ob- 
tained with model optical widths measured at 2.15 r^ 
(using a fitting function different from the Polyex model 
adopted here). This was based on the observation that 
the other RC widths examined in his work (including his- 
togram types and model widths measured at 3.2 r^) yield 
correlations with radio widths having larger zero points 
and scatters (but similar slopes). 

However, direct comparisons between different width 
estimators are not the best way to establish the most 
suitable algorithm for optical width measurement, be- 
cause they are not sensitive to the presence of systematic 
biases in the data. In this paper, we argue that model Ha 
widths measured at the optical radius (defined as a cumu- 
lative quantity and not as a fixed multiple of exponential 
disk scale lengths) provide the best match to Hi widths 
because they are the least affected by systematic effects. 
In fact, we have demonstrated that histogram widths de- 
viate with respect to measurements performed at a fixed 
radial scale in a systematic and predictable way, depend- 
ing on the shape and extent of the RCs. Next, we have 
compared model widths measured at i? pt and 2.15 r^, 
Wpe and W2.2, and shown that their difference is larger 
for galaxies with higher disk central SB. The comparison 
with Hi rotational velocities indicates that Wpe measure- 
ments are not affected by SB effects, contrary to W2.2 
ones. As shown in Figure dj the difference between Hi 
and Polyex widths plotted as a function of the i?23.5/?"d 
ratio shows no systematic trends, whereas W21 — W2.2 
increases for larger -R23.5/f*d values (i.e. higher SB ob- 
jects), varying from ^20 to ~60 km s _1 . However, even 
Polyex widths are not completely immune to systematic 
effects, as the W21 — Wpe difference (or their ratio) still 
shows a slight dependence on the radial extent of the 
observed RC, for which we provide a statistical correc- 
tion. Such dependence also exists for histogram and W2.2 
widths (see ICatinellal[200l . 

Systematic effects such as those mentioned above must 
be taken into account in order to avoid biases in the anal- 
ysis of TF surveys. As mentioned in Paper I, a width 
difference of 20 km s _1 translates into approximately a 
0.2 mag offset on the TF plane, or 500 km s _1 in peculiar 
velocity, at cz ~ 5000 km s . Thus spurious bulk flows 
of significant size can arise when heterogeneous data sets 
(e.g., optical and radio) are combined that do not share 
similar distribution properties, such as in sky projection 
or redshift. The same considerations apply to a homoge- 
neous sample of galaxies spanning a range of central SB 
and redshift, if W2.2 is used as width estimator. In fact, 
higher redshift bins would preferentially include higher 
SB objects and these require larger width corrections, as 
seen in Figure [4)3. 

4.2. Surface Brightness 

In this section we discuss in more detail the impact of 
galaxy SB and RC shape on optical width measurements. 
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Velocity widths obtained from model fits to the RCs 
depend on both RC shape and radius at which the mea- 
surement is performed. When model widths measured at 
2.15 rd and i? pt are compared, a dependence on disk cen- 
tral SB is also observed. This was shown in Figure [3] and 
is illustrated more clearly in Figure [3 where Wz.z/Wpe 
is plotted as a function of the J-band disk central SB, /j,q. 
The width ratio is approximately constant for low SB 
(LSB) galaxies 5 and decreases for high SB (HSB) ones. 
The reason for this trend is that HSB galaxies have larger 
-Ropt/^d ratios. Since the RCs are typically rising at 2.15 
ra, this implies that W2.2 reaches a smaller fraction of 
the velocity at i? op t in HSB galaxies. 

Observed rotational widths yield information on the 
total (i.e., baryonic and dark) mass contained within the 
radius where they are measured. In particular, W2.2 is 
measured at the radius where the RC of the exponen- 
tial, stellar disk would peak (i.e., where the contribution 
of the baryonic component to the total velocity field is 
maximized), whereas Wpe is measured at a larger ra- 
dius, where the contribution of the dark matter to the 
RC is expected to be larger. The dependence of the 
width ratio on SB seen in Figure [7] indicates that the 
fraction of total mass contained within 2.15 r<i decreases 
for galaxies of higher disk central SB. This yields support 
to the commonly accepted notion that in HSB galaxies 
the dark matter contribution becomes important only at 
large radii, whereas LSB galax ies are already dark matter 
dominated at small radii (e.g. JCasertano fc van Gorkoml 
119911: Ide Blok fc Bosmal2002h . In particular, the severity 
of the mass discrepancy c orrelates with central SB, b eing 
larger for lower SB (e.g.. iMcGaugh fc de Blok1ll998l and 
references therein). 

Since RC sh apes are known to depend strongly on lu- 
mino sity (e.g., iRubin et al]|1985t iPersic. Salucci fc Stell 
119961 : CGH06), it is interesting to look at the depen- 
dence of the W2.2/WPE ratio on both central SB and 
luminosity. Figure [5] shows how the ratio between the 
two model widths varies with /-band absolute magnitude, 
for galaxies of different central SB. The three SB classes 
correspond to those defined in Figure [3] (namely, the 
HSB/MSB and MSB/LSB divisions occur at /i of 17.5 
and 19.5 mag arcsec -2 , respectively). Since the shape of 
the RCs depends on luminosity and the i? pt / ?~d ratio de- 
pends on fj,Q, this plot clarifies the interplay between RC 
shape and choice of reference radius when one measures 
velocity widths from RCs. Typically, galaxies on the left 
(right) hand side of the plot have flat (rising) RCs, and 
the -RoptAd ratio increases from top (LSB) to bottom 
(HSB). For a given luminosity, the W2.2/WPE width ra- 
tio is smaller for HSB galaxies, because the RC shape is 
approximately fixed but the separation between the radii 
where the widths are measured increases. Conversely, 
along a "constant" SB line the i? op t/Vd ratio does not 
vary strongly, but the RC shape changes substantially, 
and larger width discrepancies are associated with rising 
RCs. W2.2 and Wpe show the best agreement for high 
luminosity, LSB galaxies and the largest discrepancy for 
lower luminosity, HSB objects. 

As mentioned above, the shapes of the RCs are known 

5 We use "LSB" to indicate the lower end of the SB interval 
spanned by our data. However, it should be noted that our sample 
includes very few truly low SB systems (see §[2}. 



to depend strongly on luminosity but also, as numerous 
authors have pointed out, on SB. Disentangling the rel- 
ative importance of these two quantities for the shapes 
of the observed RCs of galaxies spanning a wide range 
of properties is not a trivial task. Nonetheless, this 
would be helpful to improve our understanding of the 
systematic trends affecting optical width measurements 
discussed here. A detailed study of how RC shapes vary 
with the SB of a galaxy for the SFI++ sample is beyond 
the scope of this paper and will be explored in a future 
work. 

5. SUMMARY AND CONCLUSIONS 

In this work we have made a detailed comparison of ro- 
tational widths obtained from long-slit optical and single- 
dish 21 cm spectra. Because the combination of het- 
erogeneous data sets is necessary for the study of many 
astrophysical problems, this comparison is important in 
underscoring potential sources of bias. 

Since the choice of the best RC width for TF applica- 
tions is somewhat debated, we have considered different 
RC width estimators and discussed their limitations. In 
particular, we have studied their vulnerability to biases 
related to the shape and extent of the RCs and to the 
disk central SB of the galaxies. Our main conclusions 
are: 

(a) Optical rotational widths obtained from RC ve- 
locity histograms are affected by systematic biases that 
correlate with the shape and the extent of the RCs. The 
systematics observed in the data are well reproduced by 
our numerical simulations. 

(b) A more reliable width estimate that makes use of 
the spatial information contained in the RCs is obtained 
by fitting a function to the RC and measuring the veloc- 
ity at a fixed radial distance from the center of the disk. 
We refer to these as model widths. In this work, we used 
the Polyex model (eq. Q]) to fit the RCs. Velocity widths 
measured from the RC fit at the optical radius, i? pt, are 
referred to as Polyex widths, Wpe- 

(c) The comparison between model widths measured at 
i? opt and 2.15 rd, Wpe and W2.2, shows a dependence on 
disk central SB. Larger width differences are observed for 
galaxies of higher SB, which are characterized by larger 
Ropt/rd ratios. For a fixed SB, the W2.2/WPE ratio de- 
creases for galaxies of lower luminosity, which typically 
have rising RCs. Thus this bias occurs because RCs are 
in general still rising at 2.15 ra, and the fraction of to- 
tal mass included within that radius is smaller for HSB 
galaxies. 

(d) The best match with Hi rotational velocities is ob- 
tained with RC widths measured at the optical radius. In 
fact, W21 — Wpe is always smaller than W21 — W2.2 and, 
most importantly, is not affected by SB effects, whereas 
W21 — W2 2 is larger for galaxies of higher disk central 
SB. 

(e) When compared with 21 cm measurements, all the 
optical widths examined in this work show a dependence 
on the radial extent to which RCs are measured. For the 
model widths, larger offsets from W21 are observed for 
galaxies with Ha emission traced farther out in the disk. 
If the radial extent of the Hi emission and that of the 
detectable H11 regions are correlated (see § I3.2|) . galaxies 
with larger Ha disks and rising RCs should also have 
larger Hi extents and widths. In fact, we have shown that 
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the dependence of W21/WPE on RC extent exists only for 
rising RCs and effectively disappears for flat ones. 

Lastly, we have provided statistical corrections, param- 
eterized by the radial extent of the observed optical RC, 
that can be applied to Polyex measurements in order to 
approximate consistency with Hi rotational widths. Such 
cross-calibration has been used for the deri vation of a TF 
temp late relation for the SF1++ sample ([Masters et al.1 
2006). Failure to correct observed rotational widths for 
the systematic effects discussed here will produce biases 
in astrophysical inferences, most notably in the measure- 
ment of the peculiar velocity field and the study of the 



evolution of galactic mass-to-light ratios. These biases 
can become particularly significant when data sets are 
combined that do not share similar distribution proper- 
ties, such as in sky projection or redshift. 
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Fig. 1. — The difference between Polyex and histogram widths is plotted here as a function of RC extent. Data points are running 
averages; the bins are equally spaced, with separation 0.2 -R ma x/^?opt; the error bars are Poisson errors on the mean. Velocity widths here 
and in all other figures of this paper are expressed in units of km s — 1 . 
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Fig. 2. — The width difference in FigureQ]is plotted here for 4 classes of increasing RC slope: (a) 8 < 0; (b) < 9 < 0.5; (c) 0.5 < < 1.2; 
(d) 8 > 1.2 (where 6 is the slope of the Polyex fit between 0.5 and 1.0 i?opt; its units are km s — 1 arcsec - 1 ). The data points are running 
averages, each including 25, 100, 100, and 150 measurements in panels (a) through (d), respectively; vertical and horizontal error bars 
indicate Poisson errors on the mean and bin sizes, respectively. The data lying outside the spatial range displayed are sparse and are not 
shown. The vertical, dotted line is at Rmax=H pt, where Wpe is measured (for RCs to the left of this line, the measurement of Wpe 
requires an extrapolation of the fit). The solid lines show the results of numerical simulations, as discussed in the text. 
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Fig. 3. — The difference between model widths measured at Ropt and 2.15 is plotted here as a function of -Ropt/»*d (left panels) 
and R23.z/ r A (right) ratios. In the lower panels, the data are divided in 3 classes according to the galaxy's /-band disk central surface 
brightness: /no < 17.5 (circles), 17.5 < Ho < 19.5 (asterisks), and fj-0 > 19.5 mag arcsec -2 (triangles). In all panels, data points are bin 
averages; the errors on the mean arc not shown in (b) and (d). 
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Fig. 4. — (a) Polyex widths are compared to Hi measurements as a function of the i?23.5/»*d ratio, (b) Same as (a) for RC model widths 
measured at 2.15 rj. Data points and error bars represent running averages and errors on the mean, respectively. 
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Fig. 7. — Dependence of the W2.2/WPB RC width ratio on /-band disk central SB. Data points and error bars represent running averages 
and errors on the mean, respectively. 
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Fig. 8. — The Wb.a/WpE width ratio is plotted here as a function of the /-band absolute magnitude. The lines indicate the values of 
running averages (symbols and error bars have been omitted for clarity) for galaxies with low (dashed), medium (dot-dashed) and high 
(solid) central surface brightness, as labeled. The division in SB classes is the same adopted in Fig. [3] 
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TABLE 1 

Definitions of Radii and Velocity Widths Used in This Work 



Symbol 


Definition 




Scale length of the exponential disk from the /-band profile 




Scale length of the exponential part of the Polyex fit to the RC 


^raax 


RC extent (i.e., maximum radius up to which the folded RC is sampled) 


Ropt 


Radius encompassing 83% of the total /-band light 


^23.5 


Isophotal radius measured at an /-band surface magnitude of 23.5 mag arcsec -2 


Wg 


Difference between the 90th and 10th percentile points of the RC velocity histogram 


W PE 


Twice the velocity at R op t, measured from a Polyex fit to the RC 


W 2 .2 


Twice the velocity at 2.15 r^, measured from a Polyex fit to the RC 


w 21 


Velocity width measured at the 50% level of peak intensity on each side of the HI profile 


Note. — Throughout this work, radii 


are expressed in arcsec and velocity widths in km s -1 . The corrections applied to optical and radio widths 


are discussed in §2\ 






TABLE 2 




Width Conversion Relations 


y 


y = /(x=/? max //J op t) RCs 


W21-WPE 


= -5.21 +22.46 x -4.88 x 2 all 


w 2 i/Wpe 


= 1.004 +0.058 x -0.012 x 2 all 


W2l/Wp E 


= 0.899 +0.118 x rising 


W21/WPE 


= 1.075 -0.013 x flat 



